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Myxovirus resistance A (MxA) is an antiviral protein induced by interferon a and 6 (IFN-a, IFN-8) that can 
inhibit viral replication. The minor alleles of the —88G>T and —123C>A MxA promoter single-nucleotide 
polymorphisms (SNPs) are associated with increased promoter activity and altered response to IFN-a and 
IFN-8 treatment. Here, we demonstrate that the —123A minor allele provided stronger binding affinity to 
nuclear proteins extracted from IFN-6-untreated cells than did the wild-type allele, whereas the —88T allele 
showed preferential binding after IFN-8 stimulation. Endogenous IFN-a and IFN-§ induction can be suppressed 
in severe acute respiratory syndrome (SARS) coronavirus infection. In support of our in vitro findings, a large 
case-control genetic-association study for SARS coronavirus infection confirmed that the — 123A minor-allele 
carriers were significantly associated with lower risk of SARS coronavirus infection, whereas the —88T minor- 
allele carriers were insignificant after adjustment for confounding effects. This suggests that —123C>A plays 
a more important role in modulating basal MxA expression, thus contributing more significantly to innate 
immune response against viral infections that suppress endogenous IFN-a and IFN- induction such as SARS 


coronavirus. 


Severe acute respiratory syndrome (SARS) coronavirus 
infected 8422 people worldwide and caused 916 deaths 
[1-3]. Genetic polymorphisms have been demonstrated 
to be associated with the risk of SARS coronavirus in- 
fection. Although homozygosity of the extracellular 
neck region of the L-SIGN receptor plays a protective 
role in SARS coronavirus infection [4], it was reported 
that SARS patients who were homozygous for ICAM3 
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Gly143 (dbSNP: rs2304237) were associated with higher 
lactate dehydrogenase levels and lower total white blood 
cell counts [5]. 

Type I interferons a and 8 (IFN-a and IFN-f) can 
inhibit SARS coronavirus replication in vitro [6-10] by 
inducing antiviral proteins such as myxovirus resis- 
tance A (MxA; HGNC Symbol: MX1) [11-14]. The 
minor alleles —88T and —123A of 2 promoter single- 
nucleotide polymorphisms (SNPs), —88G>T (dbSNP: 
182071430) [15] and —123C>A (rs17000900) [16] of 
the MxA gene, have been shown with luciferase reporter 
assay to be associated with increased promoter activity 
[16, 17]. The —88G>T SNP is located in an interferon- 
stimulated response element—like sequence; the —88T 
allele makes the sequence more similar to consensus 
interferon-stimulated response element [16, 18]. With 
IFN-a and IFN-@ treatment, the reporter assays showed 
a greater increase in promoter activity for —88T but 
no significant effect for the —123C>A SNP [19]. In 
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contrast, in the absence of IFN-aw and IFN-8, —123A contrib- 
uted to increased promoter activity, whereas the difference be- 
tween —88G and —88T was small or insignificant. Consistently, 
on receipt of IFN-a treatment, the MxA messenger RNA 
(mRNA) level of healthy human peripheral blood mononuclear 
cells that have heterozygous —88GT or homozygous TT ge- 
notypes has been shown to be significantly higher than those 
that have a homozygous —88GG genotype [20]. Despite these 
studies, the binding preferences of the SNPs to nuclear pro- 
teins have not been reported. 

Although exogenous IFN-a@ and IFN-§ can inhibit SARS 
coronavirus replication in vitro [6-10], most studies report low 
induction of IFN-a@ and IFN-8 in SARS coronavirus—infected 
cells, which suggests that the virus evades the antiviral activity 
of IFN-a@ and IFN-8 [21-24]. As the —123C>A SNPs affect 
basal MxA expression without IFN-a and IFN-f induction, we 
hypothesized that this SNP is likely to be important for innate 
immune response against SARS, in which endogenous IFN-a 
and IFN-@ induction is suppressed. 

Because of their proximity (35 base pair [bp] apart), the 
—88G>T and —123C>A SNPs have been assumed to be in 
linkage disequilibrium (LD) [16, 17, 25], thus leaving the 
—123C>A SNPs less studied by comparison [16, 25]. The 
—88G>T SNPs have been studied in many diseases, including 
hepatitis B [26, 27] and hepatitis C [15, 28, 29], multiple scle- 
rosis [19, 25], and subacute sclerosing panencephalitis [17]. For 
SARS susceptibility, only the —88G>T SNPs have been eval- 
uated. Sample sizes of the studies conducted were small, and 
results were inconclusive [30, 31]. It would be important to 
clarify the contribution of both the —88G>T and —123C>A 
SNPs in relation to susceptibility of SARS coronavirus infection. 

In this study, the possible differences in binding affinity of 
these 2 MxA promoter SNPs to nuclear proteins extracted from 
IFN-stimulated and nonstimulated cells were investigated for 
the first time to our knowledge. With one of the largest col- 
lections of SARS patients and exposed control subjects, this 
genetic-association study of the —88G>T and —123C>A SNPs 
for SARS susceptibility and progression was performed on 817 
serology-confirmed SARS patients and 422 seronegative house- 
hold members of SARS patients of the Hong Kong Chinese 
population. 


MATERIALS AND METHODS 


Electrophoretic mobility shift assay. Nuclear proteins were 
extracted from the 293 cell line as described elsewhere [32] . 
Extraction of nuclear proteins from the IFN-6-stimulated 293 
cells was performed after 24 h treatment with 1000 U/mL IFN- 
Bla (R&D Systems). The probes used in the electrophoretic 
mobility shift assay (EMSA) were prepared by denaturing 2 
complementary single-stranded oligonucleotides (Figure 1A) at 
95°C for 5 min, followed by cooling at room temperature. Each 
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probe was radiolabeled by [y-”P] ATP (3000 Ci/mmol/L at 10 
mCi/mL; Perkin Elmer) by using T4 polynucleotide kinase 
(Promega), followed by purification using MicroSpin G-25 col- 
umns (GE Healthcare Life Science). 

EMSA reactions were performed using Gel Shift Binding 5X 
buffer (Promega) according to manufacturer’s protocol. In 
brief, each reaction contained 7.5 yg of 293 cell nuclear protein 
extract, 0.035 ymol/L of the labeled probe, and various con- 
centrations of unlabeled probes as appropriate and then in- 
cubated at room temperature for 30 min. Bovine serum al- 
bumin was used in place of 293 cell nuclear protein extract as 
a negative control. The reaction mixtures underwent electro- 
phoresis on an 8% polyacrylamide gel (29:1 acrylamide/bis- 
acrylamide, 1X TBE). After electrophoresis, the gel was vacu- 
um-dried. The radioactive signal was detected by exposing the 
dried gel to X-ray film at —70°C overnight. Intensity of the up- 
shifted bands was quantitatively measured with a densitometer 
[33]. The EMSAs were duplicated for the 2 SNPs in each ex- 
perimental setting. 

SARS patients and control subjects. The subjects were re- 
cruited for this study as described elsewhere [5], with approval 
from the respective institutional review boards of hospitals in- 
volved. Written informed consent was obtained from SARS 
patients who donated peripheral blood and from control sub- 
jects who donated saliva. 

A total of 817 Chinese SARS patients [5] had been recruited 
either during the outbreak in 2003 or at the SARS follow-up 
outpatient clinics in 6 hospitals in Hong Kong, namely, Pamela 
Youde Nethersole Hospital, Princess Margaret Hospital, United 
Christian Hospital, Queen Mary Hospital, Alice Ho Miu Ling 
Nethersole Hospital, and Prince of Wales Hospital. All 817 
SARS patients were confirmed to be seropositive. Their clinical 
data were retrospectively obtained from the Hospital Authority, 
Hong Kong, with permission from all attending clinicians of 
the respective hospitals. 

Household contact control subjects were members of the 
households of SARS patients and who remained unaffected and 
seronegative at the end of the outbreak. All recruited household 
contact control subjects were asked to state their relationship 
with the SARS patient(s) and other family members in the 
household. In addition to the 309 control subjects (described 
elsewhere [5]), more control samples were collected, resulting 
in a total of 422 control subjects. 

Genomic DNA from peripheral blood samples were extracted 
using conventional methods and level 3 bio-safety precautions 
in accordance with the World Health Organization and Center 
for Disease Control guidelines. The household contact control 
subjects donated saliva samples for extraction of buccal DNA 
by Oragene DNA self-collection kit (DNA Genotek). The de- 
mographic characteristic features and the clinical profile of the 
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Figure 1. Electrophoretic mobility shift assay (EMSA) for the myxovirus resistance A (MxA) —88G>T and —123C>A single-nucleotide polymorphisms 
(SNPs). A, The sequence of the probes used in EMSA. S, sense; A, anti-sense. The SNPs are underlined in each oligonucleotide for the MxA SNPs. 
B, EMSA for the —88G>T SNP, using nuclear extracts from untreated 293 cells. Comparison of lanes 3 and 4 and lanes 5 and 6 showed that the 
—88T probe bound slightly stronger than did the —88G probe. Comparison between lanes 9 and 10 and lanes 11 and 12 also showed that the —88T 
probe bound slightly stronger (also shown in the bar chart). However, the intensity of the upshifted band was low, and the competition by cold probes 
was weak. C, EMSA for the —123C>A SNP. using nuclear extracts from untreated 293 cells. The presence of 2 upshifted bands suggested that there 
might be >1 transcription factor complexes binding to the probes. The stronger band intensity of the —123A probes and the stronger competition by 
the cold —123A probe suggest that the binding affinity of —123A allele is stronger than that of —123C. D, EMSA for —88G>T SNP. using nuclear 
extracts from IFN-8—treated 293 cells. There was an upshifted band of the —88T probe, which could be competed by the cold probe. However, the 
—88G probe was not upshifted. £&, EMSA for the —123C>A SNP, using nuclear extracts from IFN-G—treated 293 cells. The upshift pattern of the 
probes is similar to that with untreated 293 cells nuclear extract, indicating that the binding affinity of the probe is not significantly affected by IFN- 
8 treatment. In the figures, the arrows indicate the upshifted bands being competed away by the cold (unlabeled) probes. The relative intensities of 
the upshifted bands were quantified, and the results were plotted as bar charts. Asterisks indicate the intensity of the upshifted band used to normalize 
that of other bands in the respective EMSA assays. 
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Table 1. Demographic Feature and Clinical Profile of Patients 
with Severe Acute Respiratory Syndrome and Control Subjects 
Who Were Successfully Genotyped 


Patients Control subjects 


Characteristics (n = 792) (n = 418) 
Age 
ean (SD) 40.31 (13.8) 43.23 (14.8) 
edian 38 45 
Range 5-88 18-83 
Sex 
ale 305 (38.5) 196 (46.5) 
Female 487 (61.5) 222 (63.1) 
ale to female ratio 253.2 2:2.3 
Treated in ICU 132 (16.7) 
Required ventilation (hypoxemic) 71 (9) 
Received steroid treatment 770 (97.2) 
Received pulse steroid/IVIG 505 (63.8) 
Death 9 (1.1) 
Number of days for hospitalization 
Mean (SD) 28.24 (18.10) 
Median (range) 23 (4-235) 
Number of days in ICU (if treated) 
Mean (SD) 2.7 (10.894) 
Median (range) 0 (0-139) 


NOTE. Data are no. (%) of participants, unless otherwise indicated. ICU, 
intensive care unit; IVIG, intravenous immunoglobulin treatment; SD, standard 
deviation. 


SARS patients and household contact control subjects who were 
successfully genotyped are summarized in Table 1. 
Genotyping. Two MxA promoter SNPs (NCBI NM_002462) 
were studied: —88G>T (rs2071430) and —123C>A (rs17000900). 
Genotyping of all subjects was performed by polymerase chain 
reaction (PCR) restriction fragment length polymorphism as 
described, except with the use of different reverse primer [28]. 
The described forward primer (5’-TGAAGACCCCCAATTA- 
CCAA—3') was used with the reverse primer designed by us 
(5'-GAAACTCACAGACCCTGTGCTGA— 3’) for PCR using 
AmpliTaq Gold DNA polymerase (Applied Biosystems) ac- 
cording to manufacturer’s instruction and the following ther- 
mal cycle conditions: initial denaturation at 95°C for 7 min, 
35 subsequent cycles at 95°C for 30 s, 58°C for 30 s, and 72°C 
for 30 s, followed by a final extension at 72°C for 7 min. The 
PCR products were digested by either 1 U of Hhal or 2 U of 
Pstl (GE Healthcare). The digested products were analyzed on 
2% agarose gel (Figure 2). The size of the undigested product 
was 296 bp. The major —88G allele creates an Hhal restric- 
tion site (GCGC), which cleaves the 296-bp fragment into 260 
bp and 36 bp, whereas the minor —88T allele is resistant to 
Hhal restriction. The major —123C allele creates a Psfl site 
(CTGCAG), which cleaves the 296-bp fragment into 225 bp 
and 71 bp, whereas the minor —123A allele abolishes the Psél 
restriction site. Ten percent of all samples were duplicated to 
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confirm the genotyping data. Control samples of known ho- 
mozygous and heterozygous genotypes and no-template control 
were included in each run of PCR and restriction fragment 
length polymorphism. The calling of genotypes on the basis of 
electrophoresis was blindly checked by 2 individual researchers 
independently. 

Genetic and risk association analyses. The genotypic and 
allelic distributions between the patient groups and the control 
groups were assessed by x’ test using SPSS for Windows soft- 
ware (version 13.0; SPSS Inc.) and logistic regression using 
STATA for Windows software (version 9.2; StataCorp LP) [34], 
where applicable. Odds ratios (ORs) and 95% confidence in- 
tervals (CIs) were used to measure the strength of association. 
Hardy-Weinberg equilibrium was assessed by a x’ test, with 
P<.05 considered to be significant. The D’ and r’ values and 
the haplotype frequencies were calculated with Haploview soft- 
ware (version 4.0; Daly Lab) [35]. 

Because a portion of our patients and control subjects came 
from the same households, logistic regression analysis with the 
robust cluster method was performed to factor in the effect of 
genetically related subjects confounding the risk association. 
The subjects were clustered according to their genetic relation- 
ship. A sample subset that contained only genetically unrelated 
subjects was also generated for analysis. Logistic regression was 
also performed to adjust for the age and sex of the subjects. 
The status (patient or control) of the subject was used as the 
dependent variable; the genotypes of either one of the SNPs, 
age, and sex were used as independent variables. 

A Bonferroni adjustment [1-(1-0.05)'"] was applied to ad- 
just the level of significance for multiple testing, the ordinary 
level of significance .05 adjusted to <.025 for n = 2, and <.0127 
for n = 4 (n, number of SNPs analyzed). 

Analysis for association with clinical outcome measures. 
Twenty-three clinical outcome measures were analyzed (Table 
2). The association of the SNPs with binomial clinical outcome 
measures was analyzed using a x” test as described elsewhere 
[5]. The numeric variables were analyzed using 1-way analysis 
of variance (SPSS software, version 13.0). 


RESULTS 


Electrophoretic mobility shift assay. EMSA was used to in- 
vestigate the binding affinity of nuclear proteins to the probes 
bearing the —88G>T and —123C>A SNPs (Figure 1A). Ex- 


This figure is available in its entirety in the online 
version of the Journal of Infectious Diseases. 


Figure 2. Electrophoresis of the digested polymerase chain reaction 
(PCR) products of PCR restriction fragment length polymorphism for the 
myxovirus resistance A (MxA) —88G>T single-nucleotide polymorphism 
(SNP) and the MxA —123C>A SNP. 
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Table 2. List of Clinical Outcome Measures 
of Severe Acute Respiratory Syndrome (SARS) 
Patients 


This table is available in its entirety in the 
online version of Journal of Infectious Diseases. 


pression of MxA in the 293 cells has been reported elsewhere 
[36] and also confirmed by reverse-transcription PCR (RT- 
PCR) in this study (Figure 3A). Using quantitative PCR 
(qPCR), we also show induction of MxA mRNA by IFN-6 
treatment (Figure 3B). To investigate the differences in re- 
sponse to IFN-6 treatment, EMSA was performed using nu- 
clear-protein extracts from the 293 cells that were either un- 
treated or treated with IFN-6 before nuclear-protein extraction. 

EMSA with nuclear extracts from the untreated 293 cells. 
The MxA —88G and —88T probes both bound with the 293 
cell nuclear proteins, resulting in upshifted bands in lanes 2 
and 8 (Figure 1B) when compared with lanes with no nuclear 
proteins (lanes 1 and 7) (Figure 1B). The binding of nuclear 
proteins to each hot probe was competed against a cold (un- 
labeled) probe of the same or the other allele in separate re- 
actions. The upshifted band of the hot —88G probe in lanes 
3 and 4 was not significantly diminished in intensity by com- 
petition with the cold —88G probe, whereas competition with 
the cold —88T probe decreased the band intensity in lanes 5— 
6 in a dosage-dependent manner. Consistently, the upshifted 
band intensity of the hot —88T probe was diminished by com- 
petition with the cold —88T probe (lanes 9-10) but not by 
competition with the cold —88G probe (lanes 11-12). Because 
a stronger band intensity indicates stronger binding, these find- 
ings suggest that the binding affinity of the nuclear proteins to 
the —88T allele was slightly stronger than that to the —88G 
allele. Nevertheless, the decrease in band intensities by com- 
petition with either cold probe was weak, indicating that the 
binding of proteins to the probes was not specific or strong. 

The MxA —123C and —123A probes also bound with the 
293 cell nuclear proteins, resulting in 2 upshifted bands (lanes 
2 and 8 in Figure 1C). The band intensities of the upshifted 
—123A probe (lane 8) was obviously stronger than that of the 
—123C probe (lane 2), indicating that the binding affinity of 
nuclear proteins to —123A was stronger. The upshifted bands 
of hot —123C were almost completely eliminated by compe- 
tition with both the cold —123C probes (lanes 3—4) and the 
cold —123A probes (lanes 5-6). On the other hand, the in- 
tensity of upshifted hot —123A probe could be diminished by 
the cold —123A probe (lanes 9-10) but much less effectively 
by the cold —123C probe (lanes 11-12). 

EMSA with nuclear extracts from the IFN-B-treated 293 
cells. Although the difference in binding affinities of nuclear 
proteins to the —88G and —88T probes is weak, the interferon- 
stimulated response element in which the —88G>T SNP resides 


suggests that the SNP may affect the binding affinities to nuclear 
proteins extracted from the IFN-@-treated 293 cells, with —88T 
binding more strongly. EMSA was thus repeated using nuclear 
proteins extracted from the IFN-6-treated 293 cells (Figure 
1D). As predicted, there is an upshifted band with strong in- 
tensity observed for the hot —88T probe (lane 6) which could 
be competed away by the cold —88T probes in a dosage-de- 
pendent manner (lanes 7 and 8). However, such strong inten- 
sity shifted band was not observed for the hot —88G probe 
(lanes 2-4). 

However, for —123C>A, EMSA that used nuclear proteins 
from the IFN-6-treated 293 cells showed a similar pattern to 
that of untreated cells; ie, the upshifted band of —123A was 
obviously stronger than that of —123C (lanes 2 and 6) (Figure 
1E), and both bands could be competed by their corresponding 
cold probes in a dosage-dependent manner (lanes 2—4 and lanes 
6-8). These findings suggest that preferential binding of nuclear 
proteins to the —88T allele was enhanced after IFN-6 stimu- 
lation, but binding to the —123C>A SNP was irrespective of 
IFN-8 treatment. 

Genotype analyses of the MxA SNPs. Our EMSA results 
suggested the MxA —88G>T and —123C>A SNPs may each 
have different effects in relation to IFN-mediated immune re- 
sponse, a known important factor in the pathogenesis of SARS. 
As in vivo evidence to support our in vitro observations, genet- 
ic association study was performed. 

Of the 1239 subjects recruited, 1210 subjects (792 patients 
and 418 control subjects) were successfully genotyped for 
—88G>T and —123C>A SNPs. Both SNPs were in Hardy- 
Weinberg equilibrium for both the patients and control subjects 
(P>.05). Of these, 989 subjects (664 patients and 325 control 
subjects) were genetically unrelated (Table 3). Logistic regres- 
sion was performed with robust cluster method on genotype 
data of all 1210 subjects to factor the impact of genetic rela- 
tionships between some subjects. Results show that both —88T- 
positive and —123A-positive genotypes were significantly as- 
sociated with decreased susceptibility to SARS coronavirus 
infection, P = .010 (OR, 0.72 [95%CI, 0.56-0.92]) and P = 
.007 (OR, 0.68 [95%CI, 0.52-0.90]), respectively, with a stron- 
ger association clearly demonstrated by the —123C>A SNP 
(Table 4). 

Analysis of the 989 genetically unrelated subjects also showed 
significantly decreased risk similar to that of SARS coronavirus 
infection, for the —88T-positive genotype (P = .004; OR, 0.68 
[95% CI, 0.52-0.88]) and the —88T allele (P = .009; OR, 0.75 


This figure is available in its entirety in the online 
version of the Journal of Infectious Diseases. 


Figure 3. Expression of Myxovirus resistance-A (MxA) and induction 
of MxA messenger RNA (mRNA) by IFN-@ in the 293 cell line. 
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Table 3. Numbers of Household Subjects (Genetically Related or Unrelated) and Independent Subjects 


Control 
No. Patients subjects Total 
Subjects with no related subject in either patient or control group, namely independent subjects 432 98 530 
Subjects belonging to households containing =2 subjects (those genetically unrelated)* 360 (232)? 320 (227) 680 (459)? 
Independent and genetically unrelated subjects 664 325 989 


NOTE. There were 267 households that contained =2 subjects. The total number of patients and control subjects who participated and were successfully 
genotyped for the myxovirus resistance A —88G>T and —123C>A SNPs in this study is 1210. 


* Numbers in parentheses indicate the genetically unrelated subjects (eg, spouses, in-laws, etc). 


[95% CI, 0.61-0.93]), as well as the —123A-positive genotype 
(P = .002; OR, 0.62 [95% CI, 0.46-0.85]) and the — 123A allele 
(P = .005; OR, 0.67 [95% CI, 0.51-0.89]), again with a stronger 
association demonstrated by the —123C>A SNP (Table 4). 
To exclude the confounding effect of age and sex, logistic 
regression analysis adjusted for age and sex was performed. The 
results confirmed the significant association remained valid for 
the —123C>A SNP, whereas the —88G>T SNP became mar- 
ginally significant (Table 5). Previous reports had suggested that 
2 SNPs of 2',5’-oligoadenylate synthetase 1 (OAS-1) gene— 


OAS-1 +590A>G (181131454) and +1295A>G SNP (1rs2660)— 
were associated with the risk of SARS coronavirus infection 
[30, 31]. Our genotyping of these 2 SNPs, however, showed no 
such association (see the Appendix, which is not available in 
the print edition of the Journal; Table 6). To explore any possible 
confounding effect the OAS-1 SNPs may confer on the risk 
association of the MxA SNPs, logistic regression adjusted for 
the 2 OAS-1 SNPs together with age and sex was performed 
(Table 7). Although the OAS-1 SNPs did not significantly affect 
the association of the —123C>A SNP (P = .002), it rendered 


Table 4. Risk Association Analyses of the Myxovirus Resistance A Promoter —88G>T and —123C>A Single- 


Nucleotide Polymorphisms 


Logistic regression* SARS Caneel x’ test 
SNP Pvalue OR (95% Cl) patients? subjects” ~~ P-value OR (95% Cl) 
—88G>T genotype 014 016° 
GG 387 (58.3) 158 (48.6) Reference 
GT 233 (35.1) 141 (43.4) 
TT 44 (6.6) 26 (8.0) 
T-carrier (GT/TT) 010 0.72 (0.56-0.92) 277 (41.7) 167 (51.4) 004° 0.68 (0.52-0.88) 
Total 664 325 
HWE P value 0.270 0.481 
—88G>T allele 
G 1007 (75.8) 457 (70.3) Reference 
T 321 (24.2) 193 (29.7) 0094 0.75 (0.61-0.93) 
Total 1328 650 
—123C>A genotype .010 .008° 
EC 533 (80.3) 233 (71.7) Reference 
CA 123 (18.5) 88 (27.1) 
AA 8 (1.2) Al (2) 
A-carrier (CA/AA) .007 0.68 (0.52-0.90) 131 (19.7) 92 (28.3) .0024 0.62 (0.46-0.85) 
Total 664 325 
HWE P value 0.764 0.173 
—123C>A allele 
C 1189 (89.5) 554 (85.2) Reference 
A 139 (10.5) 96 (14.8) .005° 0.67 (0.51-0.89) 
Total 1328 650 


NOTE. Data are no. (%) of participants, unless otherwise indicated. Cl, confidence interval; HWE, Hardy-Weinberg equilibrium; 


OR, odds ratio; SARS, severe acute respiratory syndrome. 


° Analysis of all 1210 subjects with logistic regression clustered by genetic relationship. 


z Analysis of only genetically unrelated subjects. 
° P values for overall genotypes by x” test (3x2). 
4 P values obtained against reference by x? test (2X2). 
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Table 5. Risk Association Analyses of Myxovirus Resistance A (MxA) Pro- 
moter —88G>T and —123C>A Single-Nucleotide Polymorphisms Adjusted for 
Age and Sex in Genetically Unrelated Subjects 


SNP P value OR (95% Cl) 
MxA —88G>T 
Overall genotypes GG/GT/TT comparison 044 0.78 (0.62-0.99) 
Wild type vs Fcarriers (GG vs GT/TT) 025 0.71 (0.52-0.96) 
MxA —123C>A 
Overall genotypes CC/CA/AA comparison .005 0.63 (0.46-0.87) 
Wild type vs A-carriers (CC vs CA/AA) 002 0.58 (0.41-0.82) 


NOTE. Logistic regression analysis was performed using age and sex as covariates. Only 
the P values of the genotypes of the MxA promoter SNPs are reported in this table. The P 
values of age and sex in each test were significant. Although the P values obtained for 
genotype analysis of —-123C>A SNP remained significant as in the x” test, the association 
for the overall genotype of —-88G>T SNP became marginally significant (applying a Bonferroni 
adjustment for significance; n = 2, P<.025). Cl, confidence interval; OR, odds ratio. 


the association of —88G>T SNP insignificant (P = .039), com- 
pared with the adjusted significant cutoff of .0127 (Bonferroni 
adjustment for n = 4). 

Linkage disequilibrium of MxA -—88G>T SNP and 
— 123C>A SNP. The D’ values for all subjects and the un- 
related subjects were both 1.0, whereas r* was 0.391 and 0.384, 
respectively. This indicates that the SNPs are not of perfect LD. 
The haplotype frequency analysis between patients and con- 
trol subjects (Table 8) showed that with respect to the wild- 
type C''~G* haplotype, the non—C’'”G ® haplotypes (ie, 
CT 8, AT 88 and AG *) were associated with low- 
er risk of SARS (P = .009; OR, 0.75 [95% CI, 0.61-0.93]) (Ta- 
ble 8). 

Analysis for association with clinical outcome. To check 
whether the —88G>T and —123C>A SNPs might influence 
clinical outcome, genotypes were analyzed against clinical out- 
come measures for the possible association with prognosis of 
SARS. The analyses of assisted ventilation and other outcome 
measures showed no significant association with both SNPs 
(data not shown). 


DISCUSSION 


Our EMSA data showed that, for the —88G>T SNP, preferential 
binding to the minor —88T allele was enhanced after IFN-8 
stimulation, consistent with the —88T allele being more similar 
to the consensus interferon-stimulated response element; for 
the —123C>A SNP, there was preferential binding to the minor 
—123A allele irrespective of IFN-@ treatment. These results are 
consistent with previous luciferase findings that greater increase 
in promoter activity was seen with IFN-6 treatment for the 
—88T allele but with no significant effect for the —123C>A 
SNP, whereas the —123A allele contributed to increased pro- 
moter activity in the absence of IFN-6 treatment [19]. These 
data suggest that in the absence of IFN induction, the —123C>A 
SNP may have greater effect on basal MxA levels than the 


—88G>T SNP. Induction of IFNs has been reported to be sup- 
pressed in SARS coronavirus-infected cells, suggesting that the 
virus evades the antiviral activity of IFN-a and IFN-6 [21-24]. 
Thus the —123C>A SNP may probably play a more important 
role in modifying basal MxA expression, thus contributing 
more significantly to the innate immune response against SARS 
coronavirus and in the pathogenesis of SARS. 

Age and sex are well-known confounding factors for SARS, 
because individuals aged =65 years and <18 years are associated 
with a lower risk, and the proportion of women was higher in 
SARS patients [2, 37, 38], as also reflected in our case patients 
(Table 1). When logistic regression was used to adjust for these 
confounding effects, the risk association of —123C>A SNP re- 
mained significant but the association of —88G>T SNP was 
weakened (Table 5). Thus, the association with risk of SARS 
coronavirus infection was clearly stronger for the —123C>A 
SNP. 

Despite the proximity of the —88G>T and —123C>A SNPs 
with a high D’ of 1.0, our genotype data showed their allelic 
correlation was low (r* = 0.391). The absence of A'”G * hap- 
lotype in our samples is consistent with previous reports, which 
showed the —123A allele always coexisted with —88T allele, 
but never with the —88G allele [16, 19, 25], with dominance 
of the C-'’G * haplotype. The fact that C”'” can coexist with 
T * indicates the LD of these SNPs is not complete. Given the 
differences in the contribution of the —88G>T and —123C>A 
SNPs to nuclear-protein binding (Figure 1) and MxA promoter 


Table 6. Risk Association Analyses of OAS-1 
+1295A>G and OAS-1 +590A>G Single-Nucle- 
otide Polymorphism (SNP) Using x” Test and 
Logistic Regression Adjusted for Genetic 
Relationship 


This table is available in its entirety in the 
online version of Journal of Infectious Diseases. 
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Table 7. Risk Association Analyses of Myxovirus Resis- 
tance A (MxA) Promoter —88G>T and —123C>A Single- 
Nucleotide Polymorphisms (SNPs) and the OAS-1 SNPs 
Adjusted for Age and Sex in Genetically Unrelated 
Subjects 


SNP P value OR (95%Cl) 
Overall genotype comparison 
xA —88G>T .042 0.78 (0.61-0.99) 
OAS-1 +590A>G 482 0.89 (0.65—-1.22) 
OAS-1 +1295A>G 918 1.02 (0.71-1.46) 
MxA —123C>A .004 0.61 (0.44—0.86) 
OAS-1 +590A>G 428 0.88 (0.64—-1.21) 
OAS-1 +1295A>G .850 1.04 (0.72-1.48) 
Wild-type vs minor allele- 
carrier 
MxA —88G>T (GG vs GT/ 041 0.72 (0.53-0.99) 
TM} 
OAS-1 +590A>G (AA vs .210 0.77 (0.51-1.16) 
AG/GG) 
OAS-1 +1295A>G (AA vs .862 1.04 (0.70-1.53) 
AG/GG) 
MxA —123C>A (CC vs CA/ —.002 0.57 (0.40-0.82) 
AA) 
OAS-1 +590A5G (AA vs RD 0.78 (0.51-1.17) 
AG/GG) 
OAS-1 +1295A5G (AA vs .885 1.03 (0.70-1.53) 
AG/GG) 


NOTE. Logistic regression analysis of each MxA promoter SNP 
was performed with the genotypes of OAS-1 +590A>G SNP 
(rs1131454) and +1295A>G SNP (rs2660), age and sex as covariates. 
Only the P values of the genotypes of the OAS-1 and MxA SNPs are 
reported in this table. The P values of age and sex in each test were 
significant. Although the P values obtained for genotype analysis of 
—123C>A SNP remained significant as in the x? test, the association 
for the overall genotype of —88G>T SNP became insignificant (ap- 
plying a Bonferroni adjustment for significance; n = 4, P<.013, be- 
cause 4 SNPs were being tested). 


activity (depending on whether there was prior stimulation by 
IFN-@ [19]), we have also clarified an independent contribu- 
tion of both these SNPs in association with SARS coronavi- 
rus infection. 

Previous association studies of the —88G>T MxA SNP for 
SARS coronavirus infection on the basis of small sample sizes 
[30, 31] were inconclusive. Hamano et al [30] studied 44 SARS 


patients and 103 control subjects in Vietnam and found no 
significant association, whereas the He et al [31] Chinese Han 
cohort of 66 SARS patients and 64 control subjects suggest- 
ed the —88GT heterozygous genotype was associated with in- 
creased susceptibility to SARS. Another IFN-inducible gene, 
OAS-1, was also studied on these same small population co- 
horts, and the +1295A>G [30, 31] and +590A>G SNPs [30] 
reported to be associated with SARS coronavirus infection. Our 
large-scale case-control analyses of the OAS-1 SNPs, however, 
found no significant association (Table 6). Interestingly, the 
confounding effect of these OAS-1 SNPs did not significantly 
alter the association of the —123C>A SNP, but rendered the 
association of —88G>T SNP insignificant (Table 7). 

Our genetic association analysis of probably one of the largest 
case-control cohorts demonstrates for the first time, to our 
knowledge, the significant association of the —123A-positive 
genotype with lower risk of SARS coronavirus infection (Tables 
4 and 5). The apparent association of —88T became insignif- 
icant after adjustment for the confounding effects of age, sex 
and the OAS-1 SNPs (Table 7). Because SARS coronavirus in- 
fection can suppress endogenous IFN-q and IFN-@ induction 
[21-24], the —123C>A SNP probably exerts greater effect on 
basal MxA expression in innate immune response which can 
affect susceptibility to SARS. The importance of the —123C>A 
SNP was also confirmed by consistently stronger association of 
—123C>A SNP that was not affected after adjustment of age, 
sex and genetic background of the OAS-1 SNPs. 

MxA can reduce influenza viral protein [39], and mice that 
express Mx1, the mouse homolog of MxA, have been shown 
to be protected against the highly lethal human H5N1 influenza 
virus [40, 41]. A recent report [42] showed reduced and delayed 
production of IFN-6 in H5N1 virus-infected human bronchial 
epithelial cells, suggesting that the attenuation of IFN-6 re- 
sponse in innate immune response may contribute to the vir- 
ulence of H5N1 viruses. This similarity between H5N1 influ- 
enza and SARS in the evasion of IFN-a and IFN-8 induction 
suggests that the —123C>A SNP may also be an important 
susceptibility factor for H5N1 infection, which warrants fu- 
ture study. 

In conclusion, we showed that, in keeping with previously 


Table 8. Haplotype Risk Association Analyses of Genetically Unrelated Subjects 


Haplotype (—123/—88) Patients Control subjects P value OR (95% Cl) 
CxnGae 1007 (75.8) 457 (70.3) Reference 

Od lates 182 (13.7) 97 (14.9) 

Agia. 139 (10.5) 96 (14.8) 

AAG 0 (0) 0 (0) 

Non-C'#8G-88 321 (24.2) 193 (29.7) .009° 0.75 (0.61-0.93) 


NOTE. Data are no. (%) of participants, unless otherwise indicated. Cl, confidence interval; OR, 


odds ratio. 


° P values obtained against reference by x? test (2X2). 
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reported promoter assay findings [19], the — 123A allele for the 
MxA gene provided stronger nuclear-protein binding than did 
the —123C allele in the absence of IFN stimulation, whereas 
the —88T allele was shown to be preferentially bound by nuclear 
proteins extracted from the IFN-stimulated cells. In support of 
these in vitro findings, our large genetic association study pro- 
vided in vivo data that confirm that, although both SNPs dem- 
onstrated significantly decreased susceptibility to SARS coro- 
navirus infection, the association was clearly stronger for the 
—123C>A SNP. This is in keeping with the condition of sup- 
pression of endogenous IFN induction in the innate immune 
response against SARS coronavirus infection. Knowledge of 
these differences in contribution of these 2 MxA promoter SNPs 
in the innate immune response to infection would be important 
for the assessment of one’s susceptibility to infectious diseases 
that may emerge in the future, especially those sharing a mode 
of action similar to that of SARS coronavirus. 
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